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EXECUTIVE  SUMMARY 


Space-based  radar  systems  have  been  proposed  for  future  CONUS 
defense,  surveillance,  and  tactical  battle-support  functions.  In  the 
initial  planning  of  these  systems,  however,  little  attention  was  given 
to  the  deleterious  effects  of  propagation  disturbances,  particularly 
in  nuclear  environments.  To  provide  a  means  of  evaluating  these  effects, 
a  general  propagation  model  based  on  the  SATCOM  formalism  has  been 
developed . 

This  first  Topical  report  has  been  written  primarily  to  present  the 
mathematical  details  of  the  model.  The  SATCOM  channel  model  character¬ 
izes  the  propagation  effects  in  terms  of  a  time-varying  transfer  function 
or  the  equivalent  tap  delay  line  for  a  one-way  path.  A  radar  signal 
traverses  the  disturbed  region  twice;  moreover,  the  radar  processor  must 
determine  the  angular  position  of  the  target,  which  requires  a  high-gain 
antenna  system. 

To  analyze  the  degradation  of  the  antenna  beam,  the  spatial  coher¬ 
ence  of  the  channel  must  be  specified;  moreover,  the  two-way  path  must 
be  accommodated.  In  Section  I  of  the  report,  we  develop  the  general 
form  of  a  complete  radar  channel  model.  The  propagation  effects  are 
characterized  by  a  spatially  and  temporally  varying  transfer  function. 

For  some  applications,  widely  dispersed  array  elements  may  be 
employed.  In  Section  II,  we  have  generalized  the  propagation  model  to 
accommodate  angle  variations  over  a  distributed  antenna  aperture;  how¬ 
ever,  we  have  not  yet  exercised  this  aspect  of  the  model.  Rather,  we 
have  concentrated  our  initial  efforts  on  systems  with  apertures  over 
which  variations  in  the  propagation  geometry  can  be  neglected. 

A  straightforward  application  of  reciprocity  shows  that  the  signal 
structure  of  the  radar  echo  is  the  product  of  the  disturbed  forward  and 
returned  complex  signals.  It  follows,  however,  that  the  coherence  func¬ 
tion  that  characterizes  the  average  structure  of  radar  signal  depends  on 
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the  fourth-order  complex  moment  of  the  channel  transfer  function,  which 
is  not  available  in  current  SATCOM  models. 

We  show  in  Section  III  that  the  forward  and  return  paths  are  uncor¬ 
related  under  the  strong  scatter  conditions  that  are  of  most  interest. 
This  means  that  the  complex  signal  correlation  function  of  the  radar 
echo  is  the  product  of  the  corresponding  complex  correlation  functions 
for  the  forward  and  reciprocal  paths,  which  can  be  derived  from  existing 
SATCOM  Codes.  Thus,  where  it  is  needed,  the  fourth-order  moment  can  be 
evaluated. 

To  illustrate  one  application  of  the  model,  in  Section  IV  we  compute 
the  average  beam  distortion  for  an  idealized  antenna  with  a  Gaussian  gain 
function.  We  did  this  so  that  we  could  analytically  evaluate  the  inte¬ 
grals  involved.  The  results  of  that  computation  confirm  the  following 
general  guidelines  for  use  of  the  SATCOM  model  to  analyze  radar  propaga¬ 
tion  effects: 

(1)  As  long  as  the  spatial  coherence  scale  is  large  compared  to 
the  effective  aperture  size,  beam  distortion  is  negligible. 

(2)  When  beam  distortion  is  negligible,  the  radar  signal  structure 
is  simply  the  square  of  the  signal  structure  on  a  one-way  path. 

(3)  In  general,  the  forward  and  reciprocal  paths  differ  only  by 
the  scale  factor  that  accommodates  wavefront  curvature  for  a 
finite  propagation  distance. 

These  results  are  not  strictly  new.  Indeed,  they  were  incorporated, 
at  least  implicitly,  in  early  Ballistic  Missile  Defense  radar  studies. 

The  results,  however,  have  not  been  previously  developed  within  the 
framework  of  the  SATCOM  channel  modeling  approach. 

The  continuation  of  this  effort  will  apply  the  model  to  simulate 
the  degradation  of  satellite-borne  synthetic  aperture  radars.  Because 
of  (3),  it  is  only  necessary  to  simulate  a  one-way  path.  The  complete 
radar  echo  structure  can  be  derived  by  appropriately  scaling  the  signal 
and  multiplying  the  results;  when  (2)  applies,  no  scaling  is  necessary. 
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I  INTRODUCTION 


Space-based  radar  (SBR)  systems  have  been  proposed  for  CONUS 
defense,  surveillance,  and  tactical  battle  support  functions.  In  the 
initial  planning  of  these  systems,  however,  very  little  attention  was 
given  to  the  effects  of  propagation  disturbances,  particularly  those 
generated  by  high-altitude  nuclear  detonations.  Thus,  the  objective  of 
the  research  reported  here  is  to  develop  an  accurate,  but  efficient 
propagation  model  to  evaluate  the  deleterious  effects  of  propagation 
disturbances  on  SBR  performance. 

The  approach  we  have  taken  follows  the  early  work  by  Hardin  and 
Tappart*  for  the  SAFEGUARD  system,  but  draws  heavily  on  the  more  recent 
satellite  communications  (SATCOM)  models  for  predicting  the  effects  of 
propagation  disturbances.  For  SATCOM,  a  compact  channel  model  has  been 

developed  that  characterizes  the  signal  structure  with  a  comparatively 
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small  number  of  basic  parameters.  The  analysis  of  SBR  propagation 
effects  has  some  unique  features,  however,  that  present  SATCOM  models 
do  not  accommodate. 

To  introduce  the  SBR  propagation  problem,  consider  the  simple 
geometry  shown  in  Figure  1.  A  monochromatic,  spherical  wave  emanating 
from  T  will  be  distorted  by  any  intervening  structured  ionization.  The 
resultant  distortion  can  be  characterized  by  the  complex  spatial  modula¬ 
tion  function 


ycv 


a) 


where 


C 


R1  +R2 


(2) 
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FIGURE  1  S8R  PROPAGATION  GEOMETRY 


and 


z 

P 


R1R2 
R1  +R2 


(3) 


The  spatial  variable  is  measured  perpendicular  to  the  direct  ray  and 
£  is  the  frequency.  The  subscript  p  is  used  to  denote  quantities  that 
depend  on  the  principal  raypath. 


Formallv,  h  (p  ,  z  ;  f)  is  the  random  modulation  imparted  to  a 
pi  p 

plane  wave,  as  can  be  seen  by  letting  R^  approach  infinity  in  Eqs.  (2) 

and  (3).  Two  correlation  scales,  which  can  be  obtained  from  the  SATCOM 

models,  are  important:  (1)  the  spatial  correlation,  i  ,  which  is  a 

■*  0  ■* 

measure  of  the  minimum  correlation  distance  of  h  (p,  ,  z  ;  f)  along  o  , 

pip  i 

and  (2)  the  coherence  bandwidth,  f  ,  which  is  a  measure  of  the  frequency 
-►  0 
coherence  of  h  (p  ,  z  ;  f ) . 

Pi  P 

The  channel  transfer  function  h  (?,  ,  z  ;  f)  characterizes  the 

pi  p 

propagation  effects  on  a  one-way  path.  What  remains  is  to  characterize 

the  signals  that  an  actual  SBR  system  will  encounter.  A  function 

diagram  for  a  monostatic  SBR  is  shown  in  Figure  2.  The  annenna  is 

characterized  by  its  aperture  distribution  function,  A(p^).  The  square 

of  the  spatial  Fourier  transform  of  A(p  )  gives  the  power  pattern  shown 

3  1 

schematically  in  Figure  2. 

Similarly,  the  transmitted  waveform  can  be  characterized  by  the 
complex  modulation  function,  vT(t),  or  its  Fourier  transform,  vT(f) . 
Accommodating  both  the  antenna  and  the  transmitted  waveform  gives  the 
general  representation  of  the  signal  at  R: 

vR(t)  -  ||A(pi')J’hp(Cp^-HrEt,zp;f+fc)vT(f) 

x  exp{2rift}df  exp{iAK*p^ ' }dp^'  (4) 

where  v£  is  the  apparent  transverse  motion  of  the  striations  as  seen 
from  R.  The  vector  AK  gives  the  pointing  direction  of  the  target  rela¬ 
tive  to  the  main  axis  of  the  antenna. 

The  target  will  generate  a  reflected  signal,  which,  when  received 
at  T,  admits  the  similar  representation 

VA(t)  *  VcRVV,Zp:f+fc)’R(f) 

•  exp{2irift}df  expdAK'P^'  Jdp^'  (5) 
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TARGET 


where 


C 


R 


(6) 


and  a  is  Che  scattering  cross  section.  The  apparent  velocity  of  the 
striations  is  invariant  to  a  source-receiver  interchange.  Thus,  the 
only  difference  in  structure  for  the  reciprocal  path  comes  from  the 
scale  factor,  C  ,  which  accommodates  the  finite  distance  of  the  signal 
source.  Stated  another  way,  the  angular  spectrum  is  not  invariant  to  a 
source-receiver  interchange,  as  can  be  deduced  from  Figure  1. 

The  SBR  processor  must  control  the  transmitted  waveform,  the  wave¬ 
form  repetition  rate,  and  the  antenna  to  detect  and  evaluate  targets, 
generally  in  the  presence  of  clutter  and  jamming.  Thus,  a  detailed 
characterization  of  VA(C)  essential  for  SBR  performance  evaluation 
and  the  design  of  mitigants. 

For  SATCOM  analyses  in  which  only  a  one-way  path  is  involved,  it  is 
comparatively  straightforward  to  compute  the  signal  correlation  function 


<vA(t)vJ(t+T)>  “^[(hphp'^ 


(6) 


which,  as  indicated  in  Eq.  (6),  is  a  function  of  the  two-point,  two- 

frequency  correlation  of  the  channel  transfer  function  h  .  Under 

P 

strong-scatter  conditions,  moreover,  the  signal  structure  can  be  ade¬ 
quately  modeled  by  a  Rayleigh  process.  Thus,  simulations  can  be  per¬ 
formed  without  explicit  recourse  to  the  diffraction  theory. 

It  follows  from  Eqs .  (4)  and  (5),  however,  that  for  the  SBR  problem 


<VA(t)vA(t+T)> 


(7) 


Thus,  knowledge  of  the  fourth-order  moment  of  the  channel  transfer  func¬ 
tion  is  necessary  to  evaluate  the  second-order  moment  of  vA;  moreover, 
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<£  involves  a  six -dimensional  integration  compared  to  a  three-dimensior.a 
integration  for  By  extending  our  earlier  analysis  of  the  scintilla¬ 

tion  intensity  structure  for  SATCOM,  we  were  able  to  show  that  under 
strong  scatter  conditions 


<h  h  *hRhR  >=<h  h'^hV  > 
NPPPP  PP  ^  P  P 


(8) 


so  that 


<v‘>v~>v4vr>M<%r>]  •  <5> 

Thus,  under  the  strong-scatter  conditions  that  are  of  primary  con¬ 
cern,  the  SBR  signal  structure  can  be  analyzed  by  running  the  SATCOM 
model  for  the  forward  and  reciprocal  paths  and  then  multiply  the  results 
Because  the  reciprocal  path  structure  is  simply  related  to  the  forward 
path,  however,  only  one  execution  of  the  channel  model  algorithm  is 
actually  necessary. 

If  the  coherence  bandwidth,  fQ,  is  larger  than  the  signal  bandwidth 
Eq.  (4)  simplifies  to 


vR(t)  «  vT(t) 


A(p  ')h  (Co  '  +  v  t,z  ;f  ) 
1  p  1  t  p  c 


*  exp{iAK*p^'  }dp^ ' 


(10) 


with  a  similar  simplification  for  Eq.  (5).  Similarly,  if  the  spatial 
coherence  scale,  lQ,  is  large  compared  to  the  extent  of  the  aperture 
distribution. 


vR(t)  «  vT(t)g(AK)hp(vEt,zp;fc)  (11) 

where  g(AK)  is  the  antenna  "voltage  pattern."  It  follows  that 
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(12) 


V‘> '  vr2(0  • 

which  is  Che  simplest  possible  case  of  interest. 

In  Section  II,  we  present  a  detailed  analysis  of  the  propagation 
effects  for  a  distributed  antenna  system.  With  appropriate  simplifica¬ 
tions,  Eq.  (10)  can  be  recovered  as  a  special  case.  The  purpose  of  the 
analysis  Is  to  accommodate  very  large  systems  in  which  the  propagation 
geometry  can  change  significantly  over  the  array.  In  Section  III,  we 
review  the  derivation  of  Eq.  (8).  In  Section  IV,  we  apply  the  SATCOM 
model  and  an  idealized  antenna  system  to  estimate  perturbation  levels 
at  which  Eq.  (12)  becomes  invalid. 

To  assist  the  reader,  a  notation  equivalence  table  has  been  gen¬ 
erated  for  the  natural  ionospheric  models  and  the  nuclear  effects  models. 
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II  PROPAGATION  EFFECTS  IN  A  DISTRIBUTED  ANTENNA  SYSTEM 


The  model  we  shall  consider  here  consists  of  an  arbitrary  config¬ 
uration  of  noninteracting  spherical  wave  radiators.  The  n  element 


emits  the  spherical  wave 


exp{-ikr  } 
n 


An  harmonic  time  variation  is  understood.  Let  a  rectangular  coordinate 
system  be  placed  within  the  element  cluster  with  the  z  axis  along  or 
near  the  direction  of  focus  as  shown  in  Figure  3. 

The  nC^  array  element  is  located  by  the  vector,  in  the  xyz 

system.  The  field  at  R  is  then  given  by  the  sunmation 


exp{i[4>_  -  kr  ] } 
E„(0,  $)  =“  I  A  - a - 2_ 


where 


rn  -  R[ 1  +  |5iJ  -  25in*aR(0,  $)]' 


61  - 
n  R 


If  the  phase,  is  chosen  to  cancel  rn  in  some  particular  direction, 
say  (0^,  4>p)  ,  the  array  is  focused  at  It  .  If  we  assume  that  <<  1, 
then 
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ErO, 


J>) 


exp{-ikR} 

R 


A  exp{-i[$ 
n  n 


ki  (3 
n 


5)1 


(17) 


where 


A  (0*  *)  “IlSil!2  -  61  ■ -a_(0,  $)  .  (18) 

n  /  n  n  tv 

The  quadratic  term  in  Eq.  (18)  is  often  neglected. 

If  a  striated  region  is  located  at  a  distance,  H,  from  the  array 

plane,  its  effects  can  be  formally  accommodated  in  Eq.  (5)  by  replacing 

A  by  A  h  where  h  is  a  complex  random  variable.  Because  h  is  random, 
n  n  n  n  n 

performance  measures  are  derived  by  averaging  quantities  of  interest. 

For  example,  the  power  pattern  of  the  array  can  be  computed  as 

P (0,  <P)  4  R2<  |er(0,  <$)  |  2  > 

*  H  AnAm<  h_h  *>  exp{i[($  -  kA  )  -  ($  -  kA  )  ] }  .  (19) 

*— <  nm  nm  n  n  m  m 

n,m 


To  evaluate  h  ,  consider  the  wave  propagating  along  the  ray  from 
-►  -*  n 

A Z  to  R.  Applying  the  Huygens-Fresnel  formula, 
n 


h(  P,  z  ) 
n  n  n 


ff. 

u(k  +  kTn)  exp 


f1Hn(*>  if 


*  exp{-iK»p  C  } 
r  sn  n 


die 


(2it)  ‘ 


(20) 


where  the  remaining  parameters  are  defined  as: 


r  (1>  r  (2) 

"en  (1)  +  (2) 

n  n 


(21a) 
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c 

n 


Clb) 


*■>  i  2 

H  *<“+(<  cos  >  +  <  sin  $  )“  tan  3  (21c) 

n  x  n  y  n  n 


p  ■  o  -  a  2  tan  G  .  (2 Id) 

sn  n  r_  n  n 


A  derivation  and  discussion  of  Eq.  (20)  can  be  found  in  Reference  5.  If 
Eq.  (20)  is  substituted  into  <  h^h^*  ^  and  c^e  ort!lo8onal  increments 
property  of  u(<) — the  Fourier  spectrum  of  the  wave  as  it  emerges  from 
the  disturbed  region — is  applied,  the  result  is 


<hnhm  >  *  <G(<  +  V  “*  (<  +  kTm):> 


•4‘k  ¥  -  \  ¥]} 


expS 


-i<*(p  C  -  p  C 
sn  n  sm  ■ 


•>} 


d< 


(2tt)  ‘ 


(22) 


With  straightforward  algebraic  manipulations,  it  can  be  shown  that 


C  ■  [1  -  (H/R  -  Sz  )/(cos  0  -  6z  )] 
n  n  n 


and 


where 


p  *  r  (1  -  sin  0  ) (cos  $  ,  sin  0  ) 
sn  n  n  n  n 


[cos 

♦  ] 

I  (cos 

i 

(5x  ) 

n| 

l 

>  ■  sin  0< 

(  ! 
(  “  i 

1  n 

(sin 

1  ( sin 

1 

W 

(23) 


(24) 


(25a) 
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and 


D  =  I  sin^  0  +  <5x  ^  +  5y  ^  -  261  *a„  sin  S']  *  .  (25b) 

n  l  n  ■'n  n  ^  J 

For  R  >>  X,  the  antenna  beam  width  is  such  that  51^  is  small  compared  to 
0 .  Thus , 


jcos  jcos  0| 

(sin  (sin 


(26) 


moreover,  sin  0  <<  1,  so  that 

n 


p  C  -p  C  SA£  (1-  HsecQ/R) 
sn  n  sm  m  ran 


(27) 


where 


A* 


nm 


A*  "  A*m 
n  m 


(28) 


is  the  vector  between  the  n  and  m  elements.  A  similar  analysis 
shows  that  the  propagation  terms  in  Eq.  (22)  cancel  and  that 

rf 

k_  ~  k_  ~  k  . 

Tn  Tm  T 

It  follows  that 


<h  h  *>  S 
n  m  ' 


d* 


<|u(k  +  k_) |  >  exp{-iK 'At  C}-  , 

i  nm 


-  exp{-JjD5(>(AintiC>}  ,  (29) 

where  D^(5)  is  the  phase  structure  function  for  the  ray  path  along  < 
and 
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C  -  1  -  H  sec  9/R 


(30) 


To  summarize,  the  far-field  voltage  pattern  for  a  cluster  of  non¬ 
interacting  radiators  with  an  intervening  propagation  disturbance  is 
given  by  the  expression 


Er(0,  $)  -  e—{-~lk-}  l  Anhn  exp{-i6tn-AK}  (31) 

n 

where 

AK  =■  k[a  (0,  $)  -  a_(0  ,  i  )  ]  (32) 

K  K  p  p 

If  the  reference  wave  is  omitted  and  the  summation  is-  converted  to  an 
integral,  the  integral  expression  in  Eq.  (10)  results. 

The  model  can  be  used  for  simulations  or  direct  computations.  As 
an  example,  we  computed  the  average  power  pattern  as  defined  by  Eq.  (19). 
Substituting  Eq.  (29)  into  Eq.  (19),  we  have 

P(0,  <fr)  ■  51  A  A  *  exp{-4D.  (Ai  C)}exp{-iAt  *Alt) 
n  m  o  $  nm  nm 

n,m 

*  r 

A(p^')A*(p^//)exp{-D5(ti(APiC)  } 

*  exp{-iAp^  •  AK}dp^'  dp^  ’•  (33a) 

By  changing  the  variable,  Eq.  (33a)  can  be  written  in  the  equivalent 
form 


P(0,  $)  ■  ||F(Ap’i)exp{-1sD(54 (Ao^C)  JexpC-iA^ ‘AKjdA^  (33b) 


where 
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F(APi) 


A(x  +  Ao,_/2)A(x  -  Ap^/2)dx 


(34) 


In  Section  IV  we  shall  evaluate  Eq.  (33)  for  a  simple  representative 
aperture  distribution.  The  form  of  the  structure  function  is  derived 
from  the  SATCOM  channel  model. 
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Ill  PROPAGATION  EFFECTS  OVER  A  TWO-WAY  PATH 


As  discussed  in  Section  I,  a  complete  analysis  of  SBR  propagation 
effects  must  accommodate  the  two-way  path.  As  noted  in  Section  I,  the 
only  effect  of  a  source-receiver  interchange  is  in  the  spherical  wave 
correction  factor  [cf.  Eqs.  (2)  and  (5)].  For  the  distributed  array 
analyzed  in  Section  III,  C  as  defined  by  Eq.  (30)  is  replaced  by 


CR  "  1  *  C 

■  H  sec  0/R 


(35) 


If  the  target  has  scattering  cross  section,  a,  it  follows  that  the 
reflected  signal  at  the  antenna  array  described  in  Section  II  admits  the 
representation 


v  (AK  )  -Gil  Ah  A  hRexp{i(dt  +  54  ) -AK  }  (36) 

so  *•  **  n  n  m  m  n  m  o 

n  m 

where  AK  is  the  wave  vector  for  a  target  in  the  O  ,  i  )  direction, 
o  o  o 

The  return  average  power  for  a  point  reflector  in  the  (3  ,  ;  ) 

o  o 

direction  can  be  computed  from  the  formula 


<|v  (AK  )|2>  -o£  £  AA*'A  A  *' 
a  o  f  -  n  n  m  tn 

«  «  •  tn  m  ’ 


nf nrm9mr 


x  exp{i(«4  -  Si  ,)-AK  } 
r  n  n'  o 


*  exp{i(Si  -  Si  ,) *AK  } 
tn  m  o 


«(h  hR,hh*> 
n  n'  mm' 


(37) 
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The  incident  power  and  R“  losses  have  been  normalized  to  unity. 
To  evaluate  the  propagation  term,  we  use  Sq.  (20)  whereby 


o* 

<h  h  h  h  *>  = 
n  n  m  m  x 


exp{-i[Kv  •  p  C 
n 


(l).^o(R)  *(2) 


,(R)  ^*(3) 


iW- 


P  Cyl'J  +KW'  •  p  C-Kv''p  ,C]> 
n'  m  m' 


X  exp|i[h(K(1))  -  h(ic(2))  +  h(ic(3))  -  h(K(4))] 
x  <dC(K(1))dC*(K(2))dC(K(3))d5*(K(4))> 


(38) 


To  evaluate  Eq.  (38),  we  make  the  following  change  of  variables 


2u(1) 

-  ?  r<R>  4.  t  r(R> 
*  0  C.  +  0  .c 

n 

n 

2£<2> 

-  t  c«)  -  ?  ,c<8> 

n 

n' 

2^3> 

.  t  c«>  +  ?  c<8> 

n 

n' 

2a(4) 

.  ?  c<8)  -  »  c(8> 

m  m' 


(39a) 

(39b) 

(39c) 

(39d) 


In  the  Fourier  domain,  we  have 


2|(1) 

-  K(1)  -  K<2>  *  *<3)  -  k(M 

(40a) 

2q<2) 

-  K(1>  +  *<2>  -  *<3>  -  S<‘> 

(40b) 

2*(3) 

■  *<»  -  it<23  -  k<3)  +  iE'4> 

(40c) 

2;<4) 

-  K(1)  ♦  k(2)  *  s<3>  +  *«>  . 

(40d) 
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Spatial  homogeneity  requires  that  Eq.  (38)  be  independent  of  A 

necessary  condition  is  that  the  fourth-order  moment  of  the  Fourier 
amplitudes  must  be  a  delta  function  in  the  q*'^  variable.  It  follows, 
therefore,  after  some  algebraic  manipulations  that 


(R)  h(R)*  h  h  *y 

N  n  n'  mm' 


exp{-i[;(2).;(2) 


(2)  .-t(3)  MV) 


I, -*(2)  -*(3)  .  M2)  HI)  dq '  '  dq v  ' 

exp^iq  -q  f  i,  q  ,  q  ,  q  )  — 7  2  2 

[  ’  (2t)Z  (2m)Z  (2*/ 


(41) 


In  the  spatial  domain  the  equivalent  expression  is 


<hiR)  hmV>  *  (*SH  jjfM* 


(2)  H2)"  H3)  M3)"  -(4) 


-  a 


a  -  a  ,  a 


where 


xe*p{-iS<2)".  -(3)-^|d^4)//da(3)//  (42) 


a(2)  -  (ap  ,C(R)  -  Ac  ,c)/2 
\  nn'  mm'  / 

j<3>  -(rj  ,c<R>  -rj  c) / 2 

\  nn'  nn'  • 


(ao  , 

\  nn' 


,(R) 


+  Ap 


3  ,c)/2 

mm'  / 


(43a) 

(43b) 

(43c) 


and  Apnn,  -  pn  -  pa„,  pnn,  -  pn  +  pn„  etc.  If  n  -  n'  and  m  -  m', 

a(2)  -  p(4)  -  0  and  a(3)  -  p  C(R)  -  p  C.  Finally,  if  C(R)  -  C,  Eq.  (42) 

n  m 

reduces  to  the  equation  for  the  intensity  autocorrelation  function  used 

4  -4  3,4 

in  <vino . 
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To  evaluate  M^,  we  use  the  Gaussian  phase  screen  model,  whereby 


M4(“(2)’  “(3)‘  “(4))  =  ^p[-[d(^(2)  +  ^)  -  ^ 

+  d(« 


t<2)  *  -O)  | 


3(3)  +  z(2)) 

+  D(a<2) 

.  :(1>) 

s<2>  -  ) 

+  D(a<3> 

-  J<4>)]! 

.  (44) 

For  a  power-law  phase  screen  with  a  one-dimensional  spectral  index  less 
than  2,  D(a)  =  C5^|aj2v  1  where  0.5  <  v  <  1.5.  When  all  the  appropriate 
substitutions  are  made,  it  can  be  shown  that 


<hRhR*  h  h  *> 
N  n  n  mm' 


exp  -  C 


[|ap  ,c 

L  nn' 


(R)  -?(2),2v-l 


-  a 


LJ.C2  fl ri  -  -C3)  |2v-l  ^  _  ;(3>,2v-l 


“Pi*C4*[|cS«i'  "S< 


+  Cp  .  -a 
1  n'm 


-  Cp  -  a 
1  nm 


(2) 


-  jO)|2v-l  .  |CJ  J<2)  .  ;(3)!2v-11J 
1  1  n'm'  J 


exp 


j.  i!  dj(3)dj(2) 


2  | 
P  ’ 


(45) 


where  Dp  *  p  C^  -  p  C,  etc. 
run  n  m 

2 

When  the  phase  structure  constant,  C.,,  is  sufficiently  large,  the 

-►(  2) 

inner  integral  approximates  5 (a  ),  To  prove  this  formally,  we  first 

consider  the  behavior  of  the  argument  of  the  first  exponential  as 
(3) 

a  “.  By  using  an  appropriate  series  expansion  [Eq.  (31)  in  Refer¬ 

ence  3],  we  can  show  that  the  argument  converges  to  zero  (as  long  as 

v  <  1.5)  so  that  the  exponential  is  near  unity.  For  small  values  of 
(3) 

a  ,  the  contribution  can  always  be  reduced  to  a  negligible  value  by 

increasing  C*,.  In  effect,  the  integral  over  is  equivalent  to  the 

-►(2) 

integral  over  the  complex  exponential,  which  is  6(a  ).  It  follows 

that 
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_2 
C.  -*» 

0  3 


lim<hRhR*  h  h  *>  S  exp(-  C2  J Ad  C(R)|2’’’_1] 
?  n  n7  m  m7  r[  ob '  nn7  ) 

f  _2  i  -*■  _  1 2v-l 

*  exp(  -  C.  ;  Ao  .C! 
r(  3$  mm'  ' 


(-6) 


When  Eq .  (10)  applies,  Eq.  (1)  becomes 


|va(ako) 


at  Y  A  A  *  <  >  exp{i(<$£**  -  6~l  )  *AK  }) 

l  »  n  n7  n  n7  n  n  o  / 

\nn'  / 

*(  Z  A  A  *  <  h  h  *  >  exp{i(<5i  -  61  ) -AK  }\  .  (47) 

l  ^ ,  m  m7  N  m  m7  '  r  m  m  o 
\  mm  j 


That  is,  the  power  received  can  be  computed  from  the  product  of  the 
average  transmit  and  received  antenna  patterns,  just  as  it  is  for  an 
undisturbed  path. 

The  analysis  presented  in  this  section  verifies  that  under  the  same 
conditions  the  intensity  scintillation  satisifes  the  Rayleigh  relation¬ 
ship 


<II'>  -  1  »  (<w'*>  |2  (48) 

where  I  =*  |v|  ,  the  forward  and  reciprocal  radar  paths  are  uncorrelated. 
This  has  proven  to  be  a  very  durable  relationship  for  naturally  occurring 
scintillation  with  near  unity.  Thus,  based  on  the  special  case  analyzed 
in  this  section  and  the  more  extensive  experimental  results,  it  seems  safe 
to  use  Eqs.  (8)  and  (9)  as  a  basis  for  analysis  and  simulations  of  SBR 
propagation  effects. 
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IV  APPLICATIONS  OF  THE  MODEL 


To  illustrate  the  use  of  the  model,  we  shall  evaluate  Eq.  (34)  for 
the  Gaussian  aperture  distribution  function 

A(p)  -  exp{-  p2/2r2)}  .  (49) 

o 


The  result  is 


F(p)  »  -rr2  exp{-  p2/(4r2)}  .  (50) 

o  o 

A  gaussian  aperture  cannot  be  realized  because  it  has  no  sidelobes;  how¬ 
ever,  it  can  be  used  to  estimate  the  distortion  of  the  main  beam  if  an 
appropriate  definition  of  the  aperture  cutoff,  rQ,  is  used. 

To  complete  the  model  specification,  we  must  define  the  mutual 
coherence  function  or,  equivalently,  the  phase  structure  function 
D^(£).  The  general  form  of  D^(^)  is  unwieldy;  thus,  approximations 
must  be  used.  For  example,  in  the  power-law  regime  where  the  three- 
dimensional  spectral  density  function  has  the  form 

*<K)  -  Cgq"(2  +l)  ,  (51) 

the  structure  function  admits  a  similar  power-law  representation 

v4">  •  c5V°in[2v'1'21  <52) 

2  5 

where  is  the  phase  structure  constant.  This  is  discussed  in  Rino. 

The  anisotropy  of  wave  field  is  accommodated  by  taking 
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V 


2  2 
C  Ac  -  B'Ad  Ap  +  A' Ad 
x _  x  y _ 

a'C'  -  B,2/4 


(53) 


2  6 

Wittwer  ’  has  developed  a  complete  signal  specification  for  evalu¬ 
ating  radiowave  propagation  disturbances.  Equation  (25)  in  his  1979 
report  is  equivalent  to  Eqs .  (52)  and  (53)  for  v  ■  1.5.  A  parallel 
development  for  ionospheric  applications  is  described  in  Reference  5  and 
the  references  cited  therein.  Table  1  summarizes  the  notations  used  in 
the  two  developments . 

For  data  interpretation,  we  have  found  it  convenient  to  use  the 
structure  and  phase  turbulent  strength  parameters,  Cg  and 

C  =>  r2X2!  C  G  (54) 

p  e  p  s 

where  r  (■  2.87  x  10  ^  m)  is  the  classical  electron  radius,  X  is  the 
e 

wavelength,  and  is  the  propagation  path  length.  On  the  other  hand, 
most  other  analyses  are  specified  in  terms  of  the  rms  electron  density, 

v  2, 

<AN  >  ,  and  the  rms  phase,  o...  The  interrelationships  are  given  bv 

e  oo 

the  formulas, 


C 

s 


8tt 


3/2 


2v-2 


r(v4Xj) 

r(v-i) 


(55) 


and 


r  -  2v-l  r(v44j)  2 

p  qo  r  ( v-y  ado 


(56) 


2 

Plots  of  (AN  >  versus  C  and  o..  versus  C  for  v 
N  e '  s  50  p 

Figures  4(a)  and  4(b). 


1.5  are  given  in 


25 


Table  1 


NOTATION  FOR  PROPAGATION  MODELING 


Natural 

Ionospere 

Wittwer 

Comment 

qo 

I'1 

s 

Outer  scale  wavenumber 

a 

VLs 

Axial  ratio  along  field 

b 

Lr/Ls 

Transverse  axial  ratio 

& 

9 

Orientation  of  transversal 
irregularity  axis 

G 

q  L 

0  2 

Geometric  enhancement 
factor 

Y  + 

BP 

0 

Briggs-Parkin  angle 

aqo/G 

r  t2,  2  2  .  2 

L  =■  L  (a  cos  $  +  sin  0) 
u  s 

Briggs-Parkin  factor 

9 

Zenith  angle 

9 

Magnetic  azimuth  angle 

z  sec  9 

K 

<ZT  -  zl)zl 

Z  =  - 

f  ZT 

"Reduced"  propagation 
distance 

V 

n  *  v  +  4 

Spectral  index  parameter 

*  -(2v+l) 

*AN  ~  q 
e 

.  2n 

*AN  ~  q 
e 

Form  of  three-dimensional 
spectral  density  func¬ 
tion 

The  Briggs-Parkin  angle— ’the  angle  between  the  line-of-sight  and  the 
magnetic  field — is  given  in  terms  of  V ,  9,  and  9  as 


cos  fgp  ■  cos  9  sin  9  cos  9  +  sin  9  cos  9 


26 


2 

A  good  approximation  to  is  given  by  the  formula 

,  C  , 

ci.  -  ^  (v  -  4.5v  +  5.5)  1  <  V  <  2  ,  (57) 

30  2tt  —  — 

which  is  essentially  equivalent  to  Eq.  (23)  in  Wittwer . ^  If  we  take 
v  =■  1.5  as  representative,  Eq.  (33)  is  readily  evaluated  as 


POA, 


♦a> 


4(lTr^2  <  2  2 


+  #AK2] 

y 


where 


(58) 


•js/  “ 


dS  =* 


- 


and 


i  - 

o 


(59a) 


(59b) 


(59c) 


(60) 


We  have  assumed  that  the  beam  is  focused  on  axis. 

It  is  readily  shown  that  as  -*■  »,  dS  -*■  0  and  .V  -►  #  -*■  1.  Thus, 
for  small  perturbations,  there  is  no  distortion  of  the  average  antenna 
pattern.  As  iQ  decreases,  however,  the  beam  is  broadened  and  its 
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symmetry  destroyed.  As  noted  In  Section  I,  the  main  parameter  is  the 
ratio  of  the  spatial  coherence  scale,  to  the  aperture  size. 

To  accommodate  the  temporal  structure  of  the  signal  as  well  as  bear, 
spreading,  Ac,  in  Eq.  (33b)  need  only  be  replaced  by  AOj_  -  5 1 .  Again, 
with  v  *■  1.5  the  integral  can  be  evaluated  analytically  as 

<vRvR'*>  -  p(0A»  $A)  Rjj  (6t)  exp  {-2  •n’»AK5 1 } 

x  exp{-(r  / fi.  )4(v  ,.5t)2}  (61) 

oo  err 


where  the  vector  n  is  given  by  the  expression 


( 


(r  n  r 

o  o 


(62) 


For  many  applications,  r  /«.  <1,  and  the  "aperture  smearing"  terms 

o  o 

in  Eq.  (42)  can  be  neglected.  Thus,  on  a  one-way  path  the  antenna 
pattern  and  temporal  coherence  functions  are  multiplicative;  moreover, 
the  multiplicative  property  is  not  upset  if  frequency  coherence  loss  is 
accommodated  as  well.  In  effect,  for  v  ■  1.5  the  SATCOM  channel  model 
can  easily  be  modified  to  predict  the  complex  signal  coherence  on  a  one¬ 
way  path.  For  a  two-way  path,  we  need  only  replace  C  in  Eq.  (59)  by  the 
corresponding  value  for  the  reciprocal  path  and  multiply  the  two  coher¬ 
ence  functions. 

The  analytic  form  presented  here  should  be  used  with  caution,  how¬ 
ever,  because  of  the  assumed  Gaussian  aperture  distribution  function. 

For  the  more  exacting  system  analyses,  a  simulation  based  on  the  model 
described  in  Section  II  should  be  used. 
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Code  F31 

ATTN:  TTV 

Office  of  Naval  Research 

WWMCCS  System  Engineering  Org 

ATTN 

Code  412,  W.  Condell 

ATTN:  J.  Hoff 

ATTN 

Code  414,  G.  Joiner 

DEPARTMENT  OF  THE  ARMY 

Theater  Nuclear  Warfare  Proj  Office 

ATTN 

PM-23,  D.  Smith 

BMD  Advanced  Technology  Center 

ATTN:  ATC-R,  0.  Russ 

DEPARTMENT  OF  THE  AIR  FORCE 

ATTN:  ATC-0,  W.  Davies 

ATTN:  ATC-R,  W.  Dickinson 

Aerospace  Defense  Conmand 

ATTN:  ATC-T,  M.  Capps 

ATTN 

DC,  T.  Long 

BMD  Systems  Conmand 

Air  Force  Geophysics  Laboratory 

ATTN:  BMDSC-HLE,  R.  Webb 

ATTN:  OPR,  H.  Gardiner 

2  cy  ATTN:  BMDSC-HW 

ATTN 

0PR-1 

ATTN 

LKB ,  K.  Champion 

Harry  Diamond  Laboratories 

ATTN 

CA,  A,  Stair 

ATTN:  DELHD-NW-P 

ATTN 

PHY,  J.  Buchau 

ATTN:  DELHD-NW-R,  R.  Williams 

ATTN 

R.  Babcock 

ATTN 

R.  O'Neil 

US  Army  Communications  Command 

ATTN:  CC-OPS-W 

Air  Force  Technical  Applications  Ctr 

ATTN:  CC-OPS-WR,  H.  Wilson 

ATTN 

TN 
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Air  Force  Weapons  LaCoratory 

Kaman  Sciences  Corp  | 

ATTN 

SUL 

ATTN 

T.  Stephens 

ATTN 

NTYC 

1 

ATTN 

NTN 

Kaman  Tempo 

ATTN 

B.  Gambill 

Air  Force  Wright  Aeronautical  Lab 

ATTN 

DAS  I AC 

ATTN 

A.  Johnson 

ATTN 

J.  Devore 

ATTN 

W.  Hunt 

ATTN 

W.  Knapp 

ATTN 

K.  Schwartz 

Air  University  Library 

ATTN 

W.  McNamara 

ATTN 

AUL-LSE 

! 

M. I .T.  Lincoln  Lab 

Assistant  Chief  of  Staff 

Studies  &  Analyses 

ATTN 

D.  Towle 

ATTN 

AF/SASC.  C.  Rightmeyer 

Mission  Research  Corp 

ATTN 

AF/SASC,  W.  Kraus 

ATTN 

R.  Hendrick 

ATTN 

C.  Lauer 

DEPARTMENT 

OF  ENERGY  CONTRACTOR 

ATTN 

R.  Kilb 

ATTN 

F.  Fajen 

Sandia  National  Lab 

ATTN 

R.  Bigoni 

ATTN 

D.  Dahlgren 

ATTN 

G.  McCartor 

ATTN 

Tech  Lib  3141 

ATTN 

F.  Guigliano 

ATTN 

Space  Project  Div 

ATTN 

Tech  Library 

ATTN 

D.  Thornbrough 

ATTN 

S.  Gutsche 

ATTN 

Org  1250,  W.  brown 

ATTN 

R.  Bogusch 

ATTN 

Org  4231,  T.  Wright 

Mitre  Corp 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

ATTN 

A.  Kymmel 

ATTN 

B.  Adams 

Aerospace  Corp 

ATTN 

G.  Harding 

ATTN 

V.  Josephson 

ATTN 

C.  Callahan 

ATTN 

T.  Salmi 

ATTN 

MS  J104/M,  R.  Dresp 

ATTN 

R.  Slaughter 

ATTN 

I.  Garfunkel 

Mitre  Corp 

ATTN 

J.  Straus 

ATTN:  W.  Hall  j 

ATTN 

D.  Olsen 

ATTN 

M.  Horrocks 

ATTN 

W.  Foster 

Analytical 

Systems  Engineering  Corp 

ATTN 

J.  Wheeler 

ATTN 

Radio  Sciences 

Pacify c-Sierra  Research  Corp 

BDM  Corp 

iTTN 

F.  Thomas 

ATTN 

L.  Jacobs 

ATTN 

E.  Field,  Jr 

ATTN 

T.  Neighbors 

ATTN 

H.  Brode,  Chairman  SAGE 

Berkeley  Research  Associates,  Inc 

Pennsylvania  State  University  j 

ATTN 

J.  Workman 

ATTN 

Ionospheric  Research  Lab 

ATTN 

S.  Brecht 

ATTN 

C.  Prettie 

Photometries*  Inc 

ATTN 

I.  Kofsky 

EOS  Technologies,  Inc 

ATTN 

B.  Gabbard 

Physical  Dynamics*  Inc  J 

ATTN 

E.  Fremouw 

ESL,  Inc 

ATTN 

R.  Ibaraki 

Physical  Research,  Inc  J 

ATTN 

R.  Heckman 

ATTN 

R.  Deliberis 

ATTN 

J.  Lehman 

ATTN 

E.  Tsui 

RAD  Associates  ] 

ATTN 

J.  Marshall 

ATTN 

R.  Lelevier 

ATTN 

C.  Greifinger 

General  Research  Corp 

ATTN 

R.  Turco  ! 

ATTN 

B.  Bennett 

ATTN 

H.  Ory 

ATTN 

W.  Wright 

Geo-Centers,  Inc 

ATTN 

M.  Gantsweg 

ATTN:  E.  Marram 

ATTN 

W.  Karzas 

ATTN 

F.  Gilmore 

Honeywell, 

Inc 

ATTN 

G.  Collyer,  Avionics  Dept 

RAD  Associates  1 

ATTN 

G.  Terry,  Avionics  Dept 

ATTN 

B.  Yoon  ,  j 
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SRI  International  Rand  Corp 


ATTN 

G. 

Price 

ATTN 

E.  Bedrozian 

ATTN 

R. 

Tsunoda 

ATTN 

C.  Crain 

ATTN 

0. 

Vickrey 

ATTN 

W. 

Chesnut 

Syl  vama  Systems  Group 

ATTN 

R. 

Livingston 

ATTN 

J.  Concordia 

ATTN 

D. 

Neil  son 

ATTN 

I.  Kohl  berg 

ATTN 

D. 

McDaniel s 

ATTN 

R. 

Leadabrand 

Visidyne, 

nc 

ATTN 

M. 

Baron 

ATTN 

C.  Humphrey 

ATTN 

A. 

Burns 

ATTN 

0.  Shepard 

ATTN 

G. 

Smith 

ATTN 

W.  Reidy 

ATTN 

V. 

Gonzales 

ATTN 

J.  Carpenter 

ATTN 

W. 

Jaye 

ATTN 

J. 

Petri  ekes 

4  cy  ATTN 

C. 

Rino 
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